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Correlations of ambient seismic or acoustic vibrations are now widely used to reconstruct the 
impulse response between two passive receivers as if a source was placed at one of them. This 
provides the opportunity to do imaging without a source, or passive imaging. Applications include 
terrestrial and solar seismology, underwater acoustics, and structural health monitoring, to cite 
only a few. Nevertheless, for a given set of data, correlations do not only yield the Green's function 
between the sensors. They also contain residual fluctuations that result from an imperfect time 
or source averaging that might eventually blur the images. In this article, we propose a heuristic 
model to describe the level of fluctuations of the correlations in the case of non-stationary wavefields, 
and more particularly in the case of scattering media. The work includes theoretical derivations 
and numerical simulations. The role of multiple scattering is quantitatively evaluated. The level 
of fluctuations decreases when the duration and intensity of the diffuse waves increase. The role 
of absorption is also discussed: absorption is properly retrieved by correlation, but the level of 
fluctuations is greater, thus degrading the Green's function reconstruction. Discrepancies of our 
simple model in the case of strong multiple scattering {k£* < 18) are discussed. 



I. INTRODUCTION 

Classical waves propagating in heterogeneous media 
have been subject to increasing interest during the last 
thirty years. Because diffuse waves show complex wave- 
forms due to the randomness of the medium, they have 
long been considered to be devoid of any deterministic 
information. Additionally, it is well established that 
conventional images (obtained with ultrasounds, radar, 
seismic waves...) are degraded when scattering increases. 
Nevertheless, mesoscopic physicists have demonstrated 
the existence of various wave phenomena that survive, 
and even develop, in the presence of multiple scattering 
Pj S H, S IE @1 ■ This mesoscopic approach has led to an 
incredible number of applications in optics, acoustics, 
oceanography 0] , and even seismology 0, [§| • Such appli- 
cations take advantage of multiple scattering to image, 
to communicate through, or to monitor heterogeneous 
media. 

Field-field correlation and passive imaging is a more 
recent idea that strongly benefits from the above de- 
velopments. The idea is that the correlation of fully 
diffuse wavefields recorded at two sensors yields the 
Green's function between them as if one sensor was 
a source. The connection between correlations and 
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the Green's function is not new and can be derived 
from the fiuctuation-dissipation theorer n I l0l| . But 
more recently. Weaver and Lobkis [ll|, Il2j| proposed 
another original approach that uses diffuse waves to 
reconstruct the exact impulse response between two 
sensors. Their experiments were followed by an amazing 
breakthrough in seismology p^ . fl^ , where earthquakes 
are not controlled but sensors easily handled. 

In practical applications like seismology, we mostly fo- 
cus our effort on the reconstruction of direct (ballistic) 
waves within an array of receivers. Such a reconstruc- 
tion is not trivial to obtain: correlations include the 
Green's function, plus fluctuations that are not easily 
washed out. These remnant fluctuations corresponds to 
the difference between correlations obtained after perfect 
and after imperfect averaging. These fluctuations, also 
named pseudo-noise in the following manuscript, eventu- 
ally reduce under time and source averaging. Note that 
here, the pseudo-noise is contained in the correlations 
and blurs the Green's function, this is different from the 
ambient noise constituted by natural vibrations used as 
input data for the correlations [l^, [l^ . The purpose of 
the present article is to compare the level of pseudo-noise 
to the level of the perfectly averaged correlation (the sig- 
nal). The signal and pseudo-noise terminology is chosen 
here by analogy to active source-sensor experiments. For 
simplicity, we name signal-to-noise ratio (SNR) the ratio 
between the level of signal in the correlations and the 
level of pseudo-noise. 
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In practical applications, estimating the level of 
signal and pseudo-noise in the correlations is a central 
issue. On the one hand, one has to evaluate the 
minimum amount of data needed to perform some 
passive images: how many sources to employ, what is 
the necessary record duration, what distance is best 
between receivers...? On the other hand, it would be 
a waste of time to acquire and process an excess of 
data if the signal-to-pseudo-noise ratio (SNR) in the 
correlation is satisfying for less. This SNR quantifies 
the convergence of the correlations toward the Green's 
function. Different theories have been developed in 
helioseismology [l3| and acoustics [H, [l^ to describe 
the convergence with respect to the time of integration 
T, the number of sources N, the distance between 
receivers r, the frequency / and bandwidth A/ of the 
recorded signal. These theoretical approaches assume 
stationary wave fields, and are particularly adapted 
to ambient noise records. What about non-stationary 
records, like coda waves? Theoretical and experimental 
works have demonstrated that multiple scattering plays 
a central role in the time- and space- symmetry of 
the correlations [l^, [2l|, Another point is now to 

quantify how multiple scattering affects the convergence 
of correlations. Given a set of sources and receivers, we 
will see how multiple scattering improves the estimations 
of the Green's function obtained by correlations. In 
the present article, we will propose a prediction for the 
SNR that could be used prior to an experiment. The 
formulation that we will derive applies to non-stationary 
wave fields. Our theoretical model will quantitatively 
describe how an increasing multiple scattering improves 
the convergence of the correlations. 

In section II of the present manuscript, we show an ex- 
ample of Green's function reconstruction in a multiply 
scattering medium, using finite difference simulations. 
Sections III and IV are devoted to the evaluation of the 
signal and the pseudo-noise in the correlations, theoret- 
ical predictions arc confronted to numerical simulations. 
In section V, absorption is added to our theoretical model 
for the SNR of the correlations. The last section describes 
the SNR in the diffusion approximation. 
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FIG. 1: Distribution of the N=1800 sources (dots) and the 
1120 scatterers (circles) in the 50Ao x 50Ao grid. The number 
of scatterers varies from to 1200, but the number of sources 
is kept constant, along with their positions. 



To mimic practical situations like seismology, we have 
to build an experimental configuration that provides long 
records (long lasting coda), but presents feeble scatter- 
ing attenuation between the receivers (labeled A and B) 
where passive imaging is performed. This would mean a 
very large grid with low concentrations of scatterers (say 
kt" ^ 10, with k the wave number and the transport 
mean- free path). Since large grids are very time and re- 
source consuming, we chose a configuration with a maxi- 
mum of scattering in a limited grid. To reduce the effect 
of scattering attenuation within the array of receivers, 
we also removed the scatterers from the central region, 
as can be noticed in Fig. [TJ This has the additional virtue 
of providing an easy to interpret first arrival: the direct 
wave is simply a wave propagating in a 2D homogeneous 
medium, whatever the amount of surrounding scatterers 
(waves reflected on heterogeneities arrive later). 



II. EXAMPLE OF CORRELATION AND 
GREEN'S FUNCTION RECONSTRUCTION 

First of all, let us begin with a simple illustration 
of Green's function reconstruction by correlation of 
diffuse waves. To simulate wave propagation in het- 
erogeneous open media, we have chosen to conduct 
2-D numerical experiments of acoustic waves [13] • The 
wave equation is solved by a finite difference simulation 
(centered scheme), with absorbing boundaries; the grid 
is 50 Ao X 50Ao large with a Ao/30 spatial pitch (Aq is the 
principal wavelength). 
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TABLE I: physical parameters of the simulations (units in Ao, 
Ag/To and To). 

A set of 1800 sources is randomly distributed over the 
grid and is kept fixed throughout the experiments (see 
Fig. [T]). The signal e{t) emitted by each source is a 
pulse with a center frequency /q and a gaussian enve- 



3 



Notation 


Description 


t 


time (variable) in records s{t) 


To 


central period of the records 


T 


record duration 


T 


time lag of the correlations 


Ta 


decay time of records s{t) 


Tc 


coherence time of the diffuse waves 


Ta 


absorption time 



TABLE II: Time notations in the manuscript. 



est: S = I.6A0, along with their transport cross section: 
S* = I.IAq. Table [I] summarizes the physical properties 
of the simulated media for the four numerical configu- 
rations. This includes the number of scatterers (whose 
density is n), the transport mean free path £* = 
the diffusion constant D — ^ and the Thouless time 

TD = (where is the average of the square of the 
source- receiver distance). Note that these quantities are 
evaluated under the " independent scattering approxima- 
tion" . 
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FIG. 2: Top left: source signal e{t). Top right: example of 
one waveform s{t) received in the multiple scattering medium 
(experiment number 4, kl* = 12). Bottom: intensity of the 
wavefield averaged over 1800 sources, for the three scattering 
media. 



lope (100% bandwitdh at -6dB). Waveforms s^(t) and 
s_b(<) are recorded at A and B during 200 oscillations 
Tq. The distance r between the two receivers varies from 
1 to 10 Ao. Typical waveforms e{t) and s{t) are plotted 
in Fig. [21 The long tail of the record in Fig. [2]-(b), sim- 
ilar to the seismic coda, corresponds to waves multiply 
scattered on the surrounding heterogeneities. The expo- 
nential decay of the averaged intensity (T^(i) in Fig. [21-(c) 
is clearly visible for times greater than SOTq. The decay 
time is determined by the scattering properties and the 
absorbing boundary conditions (open medium), and can 



be fitted by a(t) = ctqc 



The decay times r^r corre- 



sponding to the different configurations are reported on 
Tab. m All different time notations are also recalled on 
Tab, [m The velocity in the medium is c = Aq/q. To test 
the effect of diffusion on the correlations, we have con- 
ducted different sets of simulations with different num- 
ber of scatterers. Each set of simulation is composed of 
1800 numerical runs, one for each source (one source at a 
time). Note that this is different from correlation of am- 
bient vibriations, where sources are continously and si- 
multaneously excited. The scatterers are empty cavities 
of diameter Ao/3 randomly distributed on the grid (see 
Fig. [T]). Their scattering cross-section was numerically 
estimated in average over the frequency band of inter- 
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3: Comparison of the reference waveform (the Green's 
function) with the averaged correlation Cab{t) for r = 5Ao 
and for k£* = 12. (a) denotes the direct wave, (b) the late 
arrivals corresponding to waves scattered by surrounding het- 
erogeneities, (c) Weak fluctuations are also visible. 

The averaged correlation is controlled by three addi- 
tional independent parameters: the number of sources 
N, the record duration T and the distance r between re- 
ceivers A and B. As an example, we plot in Fig. [3] the 
correlation averaged over T=200 oscillations and 1800 
sources for r — 5Xq. This correlation is compared to 
the primitive of the impulse response obtained if A is a 
source: e(i)(8)e(i)® / GAB{t)dt [23|, where (g) for convolu- 
tion. We observe that the full waveform is reconstructed 
by correlation (a,c). Nevertheless, as the averaging is 
not perfect, weak fluctuations (b) are also visible, partic- 
ularly around the direct wave (first arrival) . The purpose 
of the following sections is to estimate the level of the re- 
constructed Green's function (the signal) and the level 
of fluctuation (the pseudo-noise). 



III. AVERAGED CORRELATION: AMPLITUDE 
OF THE SIGNAL 

A. Theory 

To start with, we consider that a source emits a broad- 
band pulse e{t) that propagates in a heterogeneous and 
scattering medium and is eventually collected at A and 
B. The record s{t) can be modeled as a non-stationary 
random signal. Its ensemble average (here, over all pos- 
sible source positions) is zero, and cr'^{t) denotes its vari- 
ance. An estimate of cr^(i) can be obtained by averaging 
over a large number of sources, as we did on Fig. O If 
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T is the duration of s{t) (the record length for instance), 
then a{t) = for t > T oi t < Q. Moreover, we assume 
that cr(i), of characteristic decay time To-, evolves slowly 
compared to the coherence time of the diffuse waves Tc 
and the central period To: 



It has been theoretically established that on average 
over the source position, the correlation of two records 
yields the Green's function between the receivers 
113, nil- If we note E {} the averaging over the source 
position we obtain: 



E{sA{tl)SB{t2)}^(JA{tl)(JB{t2)p{t2-h)® j G+git2 - h)d{t2 - h) - J G~^g{t2 ~ h)d{t2 - h) 



(1) 



r 



where G"*" and G~ stand for the causal and anti-causal 
Green's function, and p{t) the coherence of the diffuse 
waves. In the simplest approach, s{t) is modeled as a 
shot noise i.e., a series of replica of the initial pulse e(t) 
with random and independent arrival times [2^] . In that 
case, it can be shown that pir) is simply 



J e{t)e{t + T)dt 
Ie^t)dt ' 



(2) 



and its typical width Tc — J p^{t)dt is entirely deter- 
mined by the pulse shape. Nevertheless, depending on 
the scattering and absorption properties of the medium, 
p{t) might be slightly different. In particular, p{t) is 
spread when correlation of scatterers are observed in the 



J 



medium 

imA. 



This latter point is discussed in section 



For clarity, we employ in the following the notation: 
p{t)(E) [jG+{T)dT~ jG-{T)dT] = [p®G](t). Due to 
the spatial symmetry of the configuration and the ho- 
mogeneity of the source distribution, we have aA(t) = 
csit) — cr{t)- So far the statistical averages we men- 
tioned where relative to the source position, while the 
scatterers positions were fixed. If we replace the source 
averaging by an average over the scatterers positions (), 
we obtain a similar equation, except that the Green's 
function (G)^^ is now the effective medium Green's func- 
tion. Correlations would then read: 



{sA{tl)sB{t2)) = CrA{tl)aB{,t2)pit2 - h) ( 



(G)+B (<2 - ti)d{t2 -h)- I {G)^g {t2 - h)d{t2 - h) 
I 



(3) 



This latter form will not be discussed in the present 
manuscript, in the whole manuscript we do not employ 
any averaging ovr disorder, but source and/or time av- 
eraging. G~^ and G~ are therefore the exact Green's 
functions. 

In practical applications like imaging, the repartition 
of sources and scatterers is fixed. The Green's function 
between two sensors A and B is estimated by a double 
averaging : over a finite time T and over a finite amount 
of sources. This estimate writes : 



Cab{t) = T^Xl / ^A{t)sB 



{t + T)dt (4) 



a{t)<7{t + T)dt[p(S,GAB\{T)+F{T) (5) 



The ensemble average (over all possible sources posi- 
tions) of this estimate is given by the left-hand side of 



Eq. [5l it corresponds to the contribution that is useful for 
the reconstruction of the Green's function. More quanti- 
tatively, it predicts for all times r (including late and dif- 
fuse arrivals) and distances r the amplitude of the signal 
part in the correlation. The right-hand side of Eq. [5] cor- 
responds to remnant fluctuations F which are expected 
to tend to zero with increasing time T and number of 
sources N . The amplitude of this pseudo-noise will be 
evaluated in section IV. 



B. Numerical validation 

To confirm quantitatively the level of signal in the cor- 
relations (Eq.O, we plot in Fig.|4]the correlation Cab{t) 
(solid lines) obtained in the numerical simulations (aver- 
age over T = 200 oscillations and 1800 sources), and the 
theoretical expectation cr(t)a-{t + T)dt[pi^ Gab](j) 
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FIG. 4: Amplitude of the signal in the cross-correlations of 
diffuse fields. Theory from Eq. [5] (dotted line) perfectly fits 
the correlations (solid lines). Plot r — 5Ao corresponds to 
Fig. El 



(dotted lines) for increasing distances r. For simplicity, 
we zoomed into early times r where the direct wave 
is perfectly fitted. Note that reflections arriving later 
(r > 14) are also reconstructed and fitted with the 
proper amplitude, though they are not shown here (see 
the reconstruction of early and late arrival amplitude in 
Fig. [3]) . The agreement is perfect at all distances and all 
times, meaning that we have a satisfactory theoretical 
model for the signal level in the correlation. 

To summarize, we have here proposed a model to pre- 
dict the amplitude of the signal in the correlation of dif- 
fuse waves. From a practical point of view, in this paper, 
we are interested at evaluating the SNR of the correla- 
tions around the direct wave. Therefore, in the following, 
the signal level will be defined as the maximum of the 
direct wave, which in theory reads: 



Stheo{r,T) = max < / a{t)a{t + T)dt[p (g) Gab] (t) 



and in the simulated data: 

Snum {r, t) = max | Cab (t) | , 



(7) 



with the maximum taken around the first arrival: t ± 
r/c. Note that these definitions can easily extend to any 
part of the Green's function, including late reflections 
and coda waves. 



IV. FLUCTUATION OF THE CORRELATIONS: 
AMPLITUDE OF THE PSEUDO-NOISE. 



A. theory 

Fluctuations of the correlations are visible as long as 
the averaging is imperfect, and may blur the correlations 
if not reduced enough. In most applications, we seek to 
get a fluctuation level as low as possible, but in practice 
they are rarely negligible. Knowing and predicting the 
level of fluctuations will allow us to evaluate the relative 
error in arrival time for application like imaging. It will 
also allow us to interpret more clearly weak oscillations 
in the correlations that could either be reflections in the 
media (real signal) or just remnant fluctuations (pseudo- 
noise). The SNR will be deflned as the ratio between 
the amplitude of the averaged correlations, and the level 
of fluctuations. To begin, we chose to deflne the pseudo- 
noise in the correlations from their variance: 



var 



{Cab} = E {Cab(t)'} - E {CabIt)}' 



(8) 



where the bar denotes flnite source and time averaging, 
and the estimate is obtained over source position averag- 
ing. The theoretical derivation of the variance is given in 
appendix A. We assume that 1) A and B are distant by 
a few wavelengths, 2) as in section II, the decay charac- 
teristic time of a is much greater than the propagation 
time r/c, which is itself greater than the diffuse wave co- 
herence time Tc- Then, in the case of a single source, the 
variance simplifles as : 



vartheo 



a^{t)a^{t + T)dt / p'^{t)dt 



(9) 



(6) 



The theoretical SNR for one source can be deduced 
from above formulas (Eq. [S]and Eq. [S]): 
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SNRtheo{r, t) 



[p ® Gab] (t) ^ a{t)a{t + T)dt 



(10) 



r 



This equation can be generalized to TV uncorrelated 
sources, we have: 

SNRtheo{N) = SNRtheo{N = 1)Vn. (11) 

The ^/N dependency is expected even in the case of 
correlated sources, as long as the correlation is short- 
range. From this general form, it is interesting to note 
that The SNR depends directly on the Green's function 
Gab, and that it is hardly possible to predict the SNR 
without having an approximative idea of the Green's 
function. From Eqs.fTOlandfTl] we can also deduce that: 

1. The convergence of the correlation strongly de- 
pends on the envelope a{t) of the raw records 
s{t). Since diffusion strengthens late arrivals, the 
stronger diffusion the better the SNR. This is a 
central result of our paper. 

2. The SNR depends on the amplitude of the reference 
Green's function {Gab), so shorter distances AB 
and shorter times r are more easily reconstructed. 

3. The broader the spectrum (the smaller Tc), the bet- 
ter the SNR. 

4. The more sources, the better the SNR. 

All these features have been observed in previous exper- 
iments [l3, m, [13, HI] ■ Moreover, these latter equations 
can be simplified as follows. First, we assume that the 
coherence time of the scattered waves Tc — J p^{t)dt is 
simply determined by the duration in the initial signal 
e{t) (shot noise approach). Second, we assume that the 
record time T is greater than or of the order of To- (the 
characteristic decay time of a), so that a{t) takes the 
simple form of CToe"*/"^". The SNR for our simulation 
configuration now rewrites: 



Here follows a validation of our model with numerical 
simulations in the case of non-stationary diffuse wave 
fields. 



B. Numerical validation for A'^= 
variable. 
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FIG. 5: Numerical (stars) and theoretical (line) SNR versus 
the integration time T for N=l source, and r = SAq. The 
more the diffusion, the faster the convergence. Top: theory 
perfectly fits the data for ki* = 36 (left), but slightly differs 
from simulations for large record time and strong diffusion 
(tc — 0.44). Bottom: a better fit is obtained by adjusting the 
coherence time: Tc =0.44, 0.50 and 0.61, respectively. 



SNRtheo{r,T) « [p(E)Gab]{t) 



(12) 



This result holding for non-stationary wave fields com- 
pares to previous results in the case of correlation of 
ambient stationary noise (see for instance the statistical 
approach in [13, [3, [l^ and a geometrical approach in 
[23, [131): if we assume a stationary wave field (cr(t) = o-q) 
in a rather homogeneous medium, and if we approximate 
the coherence time by the inverse of the frequency con- 
tent of the source Tc oc , then we recover the previous 
prediction of the SNR for the direct wave in 2D: 



For practical reasons, we chose to define the pseudo- 
noise level as the average of the variance of the corre- 
lations before the direct arrival, and where the Green's 
function is null. This variance corresponds to the inten- 
sity of the fiuctuations around r = 0: 



Varnurn {C Ab} = VOTn 



(14) 



with the time average performed over 
— ^ + Tc<T<^ — Tc. The signal-to-noise ratio is 
eventually evaluated by the ratio between the signal 
level of the correlation, and this averaged variance: 



SNRtheo{T,r) oc 



/TA/c 
rf ' 



(13) 



SNRtheo — 



heo 



St 

^Jvartheo 



SNRn 



Snum 
yJVQjTjiiiTa 



(15) 
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As a result, SNR < 1 when fluctuations dominate, and 
SNR ^ 1 when correlations have converged to the refer- 
ence Green's function. In Fig.[5l we compare the theoret- 
ical and numerical SNR for = 1 source, for an increas- 
ing record length T, and for three different scattering me- 
dia. First of all, the SNR increases with the record length 
T. Moreover, the fit between numerical simulations and 
theoretical SNR is satisfying for k£* = 36. Nevertheless, 
for stronger diffusion {k£* < 18) and for large record 
times (much greater than the scattering mean free time) , 
the prediction for the SNR is found to be slightly inap- 
propriate: the pseudo-noise level around r = is found 
to be stronger than expected. This is clear evidence that 
our naive model of shot noise for diffuse waves can not 
account for the complexity of the field after several or- 
ders of scattering. In other words, we suspect that there 
exists some remnant time-correlations in s{t) that appear 
when scattering is increased. This phenomenon was ob- 
served in time-reversal experiments where a saturation 
of the SNR was observed when scattering was strongly 
increased [3l[. In our case, we imagine two possible in- 
terpretations for this phenomenon. First, the scattered 
wavefield can excite the same scatterer placed perpen- 
dicularly with respect to A and B several times, which 
results in coherent correlations around t = 0. As the 
scattering medium does not move (no averaging over dis- 
order), this contribution hardly vanishes under time av- 
eraging. Second, there exist recurrent scattering (closed 
loops) in the medium that can be excited several times by 
the same source. Both interpretations mean that s{t) can 
not be modeled as a shot noise with a constant coherent 
time. Correlations between the arrival times (scattering 
paths) induce an increase of duration of p(t), which also 
increases Tc- Fits of the numerical result yield the fol- 
lowing values : for k£* ~ 18 the best fit is obtained for 
Tc = 0.50, and for kt = 12: = 0.61. 

Moreover, we clearly see from Fig. [5] that the SNR for 
one source and 200 oscillations is lower or of the order of 
one, which means that correlations have not (yet) con- 
verged: in these conditions the direct wave between A 
and B for r ~ 5Ao is hardly reconstructed. Additional 
averaging over sources is needed. This is addressed in the 
next subsection. 



C. Numerical validation for T=200 oscillations, N 
variable. 




FIG. 6: Two possible configurations of paths that result in 
time-correlation in records s{t). Solid and dashed arrows 
stand for waves at difi'erent time, (a) The same scatterer 
is excited several times, (b) Existence of recurrent scattering 
(double arrow and black scatterers). 
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FIG. 7: SNR for r = 5Ao estimated for an increasing num- 
ber N of sources and for increasing scattering in the medium 
for various k£* ranging from oo (homogeneous medium) to 
13 (strong multiple scattering). The stronger the diffusion, 
the better the correlations. For A'^ > 100, a discrepancy is 
observed and might be due to long-range correlation in the 
disordered medium. 

We now evaluate the effect of increasing the number 
of sources on the SNR. In order to have comparable 
records (with similar statistical content), 1) sources 
are placed at random (the position of the sources are 
spatially uncorrelated); 2) each couple of record SA,B{t) 
originating from the same source is normalized by the 

maximum of | ^ J s\{t)dt, ^ s\{t)dt^ . As diffuse 

waves have a coherence length of the order of half 
the wavelength, assumption 1) is likely to be slightly 
inappropriate here. The residual spatial correlation 
between the sources' positions has been evaluated in 
appendix [Cl and its effect is found to be small for the 
1800 sources of our numerical experiments. 

We plot in Fig.[7]the SNR versus the amount of sources 
used in the averaging, for four different experimental con- 
figurations (see Tab.[l|. After 1800 sources, we obtain a 
satisfying SNR (> 10) in all cases. Up to 100 sources, 
the SNR grows like a/TV, which means that records can 
be considered as uncorrelated. As for one source, the 
SNR is always better when scattering is stronger (shorter 
k£*). It is particularly interesting to note that all scat- 
tering media provide better SNR than the one obtained 
in the homogeneous one (without scatterers). However, 
for more than 100 sources, the SNR is no longer a power 
law. The rate of convergence with the number of sources 
is slower than expected, which means that the contribu- 
tions from different sources are not totally independent. 
As a simple picture, we can again invoke that long range 
correlations might exist in the scattering media. As a 
result, the SNR in homogeneous medium eventually goes 
beyond the SNR with diffuse waves for more than 400 
to 900 sources, depending on the amount of scatterers. 
This is a clear indication that the fluctuation-dissipation 
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theorem does not totally apply to the reconstruction of 
the Green's function by correlation of diffuse waves. In 
other words, time averaging, source averaging or ensem- 
ble averaging can not be simply interchanged in the case 
of strong multiple scattering. 



V. ROLE OF ABSORPTION 




z -1 

20 40 60 80 100 




20 40 60 80 100 

time T (|J. s) 



FIG. 8: Solid line: correlation without (top) and with (bot- 
tom) absorption. Dotted line: theoretical decay. The upper 
plot is equivalent Fig. |3l Both amplitude and phase are re- 
constructed, even when absorption is present. 
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FIG. 9: SNRa to SNR ratio, evaluated for r = 5Ao, kl* = 36 
and Absorption time here is Ta = 40 

In practical applications, the medium is often (at least 
slightly) absorbing. Several authors have worked out the 
effect of absorption in the correlationjU [2^, [s^l- They 
proved that in lossy environment, the Green's function 
(including absorption) is still retrieved by correlation. 
Nevertheless the role of absorption on the convergence 
of correlations is not established. We propose in the 
following to adapt the theory presented in previous 
sections to the case of weakly absorbing media. 

First of all, let us consider that the Green's function 
with absorption is simply connected to the lossless 
Green's function G as Gj^g{T) = GAB{T)e~'^^'^'' with 
the absorption time. In our simple model, absorption 



does not depend on the frequency. The variance of 
the field is also affected: aa{t) — cr(i)e^*/'^" . First we 
want to check that the Green's function is actually 
retrieved by correlation in the numerical simulations. 
To that end, the same numerical simulations presented 
above are slightly modified to account for the absorp- 
tion. Data are processed again, an example of averaged 
correlations is presented in Fig. [S] As expected, the 
reference Green's function [p (g) G°] (r) is retrieved both 
in phase and amplitude, meaning that the absorption is 
actually reconstructed in the correlation. Nevertheless, 
fluctuations are now stronger than in a medium free of 
absorption. 

The effect of absorption on the level of fluctuations 
and the SNR is quantified now. Theory from Eq. [TUKJTT] 
can be adapted to take absorption into account. For 

simplicity, we assume that aa{t) = croe^*(~^~), and 
again T ^ Ta and Tq . Then 

(16) 

This formulation contains two contribution. First, ab- 
sorption reduces the effective length of the coda record. 
Second, the direct wave is attenuated by a factor e~"^. 
As a result, the SNR decreases with increasing absorp- 
tion. In other words, the stronger the absorption the 
lower the SNR. In Fig. El we plot the ratio SNRa/ SNR 
around the direct wave versus the number of sources 
iV, for the experimental configuration number two (see 
Tab. IJ). For the physical values retained here (r — 5Ao, 
To- = 110 and Ta — 40), theory from Eq. [11] predicts a 
ratio of SNR of 45%, which is of the order of experimen- 
tal results (the average value is 60%). The discrepancy 
comes from the fact that the actual form of a is not a 
simple exponential decay. 



VI. EXTRAPOLATION TO FULLY DIFFUSIVE 
MEDIA 

All the above numerical results have been obtained in 
a configuration where the direct wave within the array 
of receivers does not suffer from scattering attenuation: 
no scatterers were placed in the central area of the sim- 
ulation. In real experiments, scatterers are likely to be 
distributed everywhere in the medium. We propose in 
the following to extend our model to this latter configu- 
ration, under the diffusion approximation. In this case, 
the average envelope of the wave field received at a dis- 
tance R in an infinite medium is: 

(j(R,t) = \ — — e~iDt . 
^ ' ^ y ADt 

This model is particularly valid when sources are far away 
from the receivers, meaning R ^ £. We additionally 
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assume that the distance r between A and B is much 
smaller than R so that the envelope is the same for A 
and B. The effective time length of one record is similar 
to the Thouless time: 



TD 



AD' 



therefore we would expect the SNR to grow like: 



SNR(xGABir,T) 



Ntd 



On the one hand, scattering increases the Thouless 
time Tu, thus increasing the SNR. On the other hand, 
the direct wave is attenuated by scattering. These two 
effects are in competition. For the direct wave, we can 
approximate the attenuation by the product of the ge- 
ometrical spreading and the diffusion attenuation. For 
simplicity, we also assume that £ — £*. The correspond- 
ing SNR then writes: 



SNR(x 



which formally means that the SNR increases with in- 
creasing scattering ii r < £* . Again, this formal result is 
a central point of our article. Note that according to this 
equation, the SNR saturates, then decreases for strong 
scattering or large distances (r > £*). This is due to 
the scattering attenuation of the reconstructed Green's 
function (in the correlation). The best SNR is therefor 
expected for r = £* . 



VII. DISCUSSIONS AND AMPLITUDE 
RECONSTRUCTION 

In the present article, we correlated diffuse waves to re- 
construct the Green's function between passive sensors. 
This processing was carried out without any transfor- 
mation of the raw record s{t). The only modification 
that we realized was normalising each couple {sa,s_b} 
by the maximum of { / s\{t)dt, J s^(i)d<}. In particular, 
we did not change the frequency content of the records 
(no whitening), nor change the amplitudes of particu- 
lar waveforms within s(t). This has the virtue of recon- 
structing the Green's function both in phase and ampli- 
tude. Practically speaking, this allowed us to retrieve 
the geometrical spreading of the wave, along with the 
absorption, and confirms the possibility of mapping the 
attenuation of the medium under investigation "without 
a source". Nevertheless, in most practical applications, 
the amplitude of the raw records s{t) is strongly modi- 
fied. Several non-linear transformations have been pro- 
posed in the literature: 1-bit processing (retain only the 
sign) , adaptive gain constant (dynamic renormaliza- 
tion that compensate the coda decay), clipping [sS], to 



cite only a few. In our model, this would be more or less 
equivalent to compensate the decay of cr, thus increasing 
the SNR. On the one hand, these latter processing were 
shown to greatly improve the final tomographic images: 
they are particularly adapted when the phase (arrival 
time) of the wave transports the quantity of interest, like 
in tomography. On the other hand, these processing are 
likely to degrade the reconstruction of the amplitude, as 
it was observed in more recent experiments [SJ, . We 
therefore strongly suggest that pre-processing treatments 
like 1-bit or whitening be used only when the reconstruc- 
tion of the phase is targeted. 



VIII. CONCLUSION 

To conclude this article, we have reported on the con- 
vergence of correlations toward the Green's function in 
the case of independent sources in a multiple scattering 
environment. In section II, we have presented numerical 
simulations and retrieved the Green's function (in phase 
and amplitude) between passive sensors. As a simple 
model, we chose to describe the coda as a superposition 
of waves arriving at random time. This offers the op- 
portunity to develop, in sections III and IV, a simple 
estimation for the mean and the variance of the correla- 
tions, from which a SNR is derived. This SNR quantifies 
the convergence of the correlations toward the Green's 

function. The SNR was found to increase like ^Z^^, 

where N is the amount of sources, td is the effective du- 
ration of the coda (the Thouless time), and Tc represents 
the duration of the source. A central point of our paper 
is that multiple scattering (large Tq- and td) provides a 
better SNR in the correlations. This results in a fascinat- 
ing observation: instead of blurring the images as is the 
case in conventional techniques, disorder here improves 
the quality of passive images. On the contrary, absorp- 
tion was found in section V to reduce the SNR and, as 
always, is a limiting factor. Our theoretical model for the 
SNR was confronted to finite difference numerical simula- 
tions. This model was found to be valid when scattering 
is limited {k£* > 18). Nevertheless, when scattering is in- 
creased, our naive description of coda waves was found to 
be slightly inappropriate: we point out short- and long- 
range correlations of diffuse wave paths as a probable 
candidate to explain the discrepancy between our theo- 
retical model and numerical simulations. 
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APPENDIX A: VARIANCE OF THE 
CORRELATIONS FOR ONE SOURCE 

To evaluate the level of fluctuations in the correla- 
tions, we evaluate the variance of the correlations for one 

I 



This central formula is similar to previous results by 
Derode et al in the case of TR and by Sabra et 
al. in the case of stationary noise correlations [ssj. Nev- 
ertheless, we can go beyond their work after a series 
of additional reasonable assumptions. 1) In a dilute 
medium {k£ » 1), [p (8) G^a] (<?) and [p Gbb] (q) are 
constituted by the source autocorrelation (a peak around 
q = 0) followed by rapidly decaying contributions (re- 
flections from surrounding heterogeneities). These con- 
tributions can be neglected: G{t) « which implies 
[p Gaa] (q) [p 8) Gbb] (q) « p'^iq)- 2) The variance a of 
the record evolves with a characteristic time of To- much 
greater than the diffuse wave coherence time Tc {tc ^ T^r), 

I 



source, averaged over the source position. Let's start 
with a simple calculation, assuming that s{t) is gaussian 
and using the moment theorem: 



(Al) 

(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



so that a{e + q) ~ a{e) and a{0 + q + T)fv a{0 + t). 3) 
E A — B ^ X, then \Gab\ |Gaa|- These assumptions 
greatly simplify the above equation that now reads: 



var[CAB{T)]- J a^i0)<J^{e + T)de J p\q)dq 

APPENDIX B: VARIANCE OF THE 
CORRELATIONS FOR N SOURCES 

We now take into account an averaging over N sources, 
each source is labeled i or j: 



E^CABirf} = J J E{sA{ei)sB{ei + T)sA{92)sB{e2 + T)}deide2 
= jj E {sA{ei)sB{ei + r)} dOiE {sA{e2)sB{e2 + r)} de2 
+ JJ E{sA{ei)sA{e2)}E{sBi0i + T).SBi02 + T)}de2dei 

+ JJ E{sA{ei)sB{e2+T)}E{sAi92)sB{0i+r)}de2dei 

The first term is the intensity of the correlations. Once eliminated, only the variance remains: 

var {Cab{t)^ = JJ a{ei)a{e2)a{ei + r)a((?2 +t)[p® Gaa] (Oi - 62) [p ® Gbb] (Oi - 02)d^?2d^i 

+ JJ <j{0i)a{ei + T)a{e2)a{e2 + r) [p Gab] {Bi - ^2 + r) [p ® Gab] {02 - Oi + T)de2d0i. 
We replace q by ^2 — ^^i : 

var |Gas(t)} = JJ cr{6i)ai9i + q)cr{9i + t)<7{9i +T + q)[p® Gaa] (q) [p ® Gbb] iq)dqd6i 

+ cr(6'i)f7(e'i +T + q)a{ei + q)cr{9i + r) [p (g) Gab] (r + q)[p(S Gab] (t - q)dqdei . 



var[CAB{r)} // ^ {■'^A(^?i)sii(^2)} {^(^i + ^)4(^2 + r)} dM^i (Bl) 

+ E[s\{9,)s'B{(^2 + r)]E[s\{e2)s^B{0i+r)]d92de^. (B2) 

I 

These summations can be split into two contributions T)d9 J p'^{q)dq, and directly derives from to the above 
i = j and i ^ j. The first one yields cr'^{9)a'^{6 + 
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case N=l. The term i ^ j exhibits cross-correlations 
^{■^^■5/1} ^^"^ ^{'^A'^'b}' which are neglected as long 
as sources are dilute enough. This term i ^ j con- 
tains short- and long-range correlations. The short- 
range correlation in diffuse media takes the usual form of 
sinc{kdij)e~'^^^ , with the distance between source 
i and j. This correlation ranges over one wavelength. 
Its role is developed in the next section of the appendix 
and can not account for the whole discrepancy between 
observed and theoretical SNR. When diffusion increases, 
their also exists some long-range contributions. These 
latter phenomena will be subject to further investiga- 
tions. 



APPENDIX C: SPATIAL CORRELATION OF 
SOURCES 

We evaluate here the number of uncorrelated sources 
when N sources are chosen at random but have are spa- 



tially correlated over A/2, which is typically the case 
inside multiple scattering media (short-range correla- 
tion) [23|. We assume that the choice of N locations 
is random over M sites. The size of the site is deduced 
from the coherence length of a diffuse field, which result 
in a coherence area of a source of the form 7rA^/4. In the 
case of a 50A x 50A large grid, we get M « 3183 uncor- 
related sites. The amount of uncorrelated sources A^e// 
is therefore lower than N, as two sources can be chosen 
at the same site. We start with Neff{l) = 1. Then, we 
iterate: 

For the 3183 independent sites, we get N^jf = 1377 
independent sources, which result in lowering the SNR 
of 23%. Short-range correlations can therfore not solely 
explain the discrepancies between the theoretical SNR 
and the SNR form numerical simulations. 
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